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^ ' Abstract 

c ' 

Context. Determining the mass of stars is crucial both to improving stellar evolution theory and to characterising ex- 
oplanetary systems. Asteroseismology offers a promising way to estimate stellar mean density When combined with 
accurate radii determinations, such as is expected from GAIA, this yields accurate stellar masses. The main difficulty 
is finding the best way to extract the mean density of a star from a set of observed frequencies. 

Aims. In this paper, we seek to establish a new method for estimating stellar mean density, which combines the sim- 
plicity of a scaling law while providing the accuracy of an inversion technique. 

Methods. We provide a framework in which to construct and evaluate kernel-based linear inversions which yield directly 
the mean density of a star. We then describe three different inversion techniques (SOLA and two scaling laws) and 
i-C ■ apply them to the sun, several test cases and three stars, a Cen B, HD 49933 and HD 49385, two of which are observed 

by CoRoT. 

Results. The SOLA (subtractive optimally localised averages) approach and the scaling law based on the surface correct- 
ing technique described by Kjeldsen et al. (2008, ApJ, 683, L175) yield comparable results which can reach an accuracy 
of 0.5 % and are better than scaling the large frequency separation. The reason for this is that the averaging kernels 
from the two first methods are comparable in quality and are better than what is obtained with the large frequency 
separation. It is also shown that scaling the large frequency separation is more sensitive to near-surface effects, but is 
much less affected by an incorrect mode identification. As a result, one can identify pulsation modes by looking for an I 
and n assignment which provides the best agreement between the results from the large frequency separation and those 
from one of the two other methods. Non-linear effects are also discussed as is the effects of mixed modes. In particular, 
it is shown that mixed modes bring little improvement to the mean density estimates as a result of their poorly adapted 
kernels. 
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' 1. Introduction sitive indicator of the mean density of a s tar. Current space 



X 



missions CoRoT (e.g . iMichel et all 120081 ) and Kepler (e.g 
Determining accurate stellar mass is crucial to several do- |Kieldsen et ail [2011 are measuring stellar pulsation {re- 
mains in astrophysics. Indeed, mass plays a dominant role quenc i es to high accuracy. When combined with an inde- 
in the evolution and final fate of stars. It is also a key p endent determination of the radius, this yields the stel- 



parameter when characterising exoplanetary systems since lar mass Than ks to very precise parallax measurements, 

53 ■ the masses of exoplanets, essential for determining whether radii accurate to 2 % are expected from the forth coming as- 

they are rocky or gaseous in nature, are in general deter- trometric mission GAIA (iPerrvman et al.ll200l"ri . The main 

mined with respect to that of the star. However, in spite challenge is then finding the best way of extracting the 

of its importance, stellar mass is usually difficult to obtain mean density of a star from a given set of frequencies, 
accurately for single stars and is frequently model depen- 
dent. Various methods have been devised or proposed for 

A method which has often been used for estimat- estimating stellar mean density from pulsation frequen- 

ing stellar mass is to compare the position of a star in cies. These include using sim ple scaling la ws based on 

an Hertzsprung-Russell (HR) diagram with evolutionary the large frequency se paration (lUlrichl fl 986) or other fre- 



tracks for stellar models of different masses. The main dim- quency combinations (jKieldsen et al.l 12008 ). searching for 



culties with this method is its model dependence, the large the best fitting model within a grid of models or a param- 
error bars, and the fact that there are regions in the HR di- eter space and retaining its mean density (e.g. Bazot et al.l 



agram in which more than one stellar evolution track goes I2008L iMetcalfe et alJl201(iL iKallinger et al.l 120 ldh . and per 
through each and every point. A promising alternative is forming a full structural inversion to determine the den- 
to use asteroseismology. Indeed, acoustic modes are a sen- sity profile, then integrating it to obtain the mean density. 
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The structural inversion may be carried out using linear 
or non- linear techniques (e.g. [Goueh fc Thompsonl Il99ll 



lAntialll996l iRoxburgh fe Vorontsovll2002D . The scaling re- 
lations have the obvious advantage of simplicity over the 
other methods and require fewer frequencies than a full 
structural inversions, but are less accurate. Although more 
expensive computationally, searching for best fitting mod- 
els in a grid is also straightforward and usually more pre- 
cise, but is model dependent, so that the conclusions are 
only as firm as the input physics. Full structural inversions 
have the advantage of being less model dependent, but usu- 
ally involve a number of free parameters which need care- 
ful adjustment and can fail to converge when non-linear. 
However, with the current missions CoRoT and Kepler and 
future asteroseismic projects, there is a growing need to pro- 
vide accurate seismic stellar masses for a large number of 
stars in an automatic way. An ideal method would combine 
the advantages of the above techniques while avoiding their 
weaknesses. 

In this paper, we focus on kernel-based linear inversion 
methods which yield directly the mean density without in- 
verting for the full density profile. As such, this represents 
a first attempt at constructing methods which combine the 
simplicity of scaling laws with the accuracy of full struc- 
tural inversions. Although the results are not as accurate 
as initially hoped for, this approach gives insight into the 
precise relation between the mean density of a star and 
pulsation modes, and provides a better way of evaluating 
scaling laws. Section [2] establishes a general framework in 
which to construct and evaluate linear inversion techniques 
which yield the mean density directly. The following section 
then describes various linear inversion techniques, including 
linearised scaling laws. This is followed by a section which 
shows how to extend these techniques beyond the linear 
regime. Sections [5l |5] and [7] applies these different meth- 
ods to the sun, to test cases using a grid of models, and 
to three stars, two of which were observed by CoRoT (HD 
49933 and HD 49385). Finally a discussion concludes the 
paper. 

2. General aspects 

2.1. Starting point 

As for most inversion methods, the methods described be- 
low start with an observed star and a reference stellar 
model. At this point we will assume that the star and the 
model have the same radius, but will relax this constraint 
later on. The structure of the reference model needs to be 
sufficiently close to the structure of the star so that the 
link between the relative frequency differences, Sv n £/u n i, 
and the differences in the stellar structure can be described 
by a linear relation: 



f K^x^dx- 
Jo P 



r 1 K f up (x) s -^dx 

Jo 1 1 



(1) 



where p is the density, Ti the first adiabatic exponent, 
x = r/R the fractional radius, n the radial order, £ the 

-Fsurf Oral) an g^hoQ 



harmonic degree, R the radius and 



surfa ce term (e.g. IChristerise n-Dalsgaard fe Berthomieul 
Il99lh . The function F surf is assumed to be a slowly varying 



funct ion of v only, an d Q n i represents the mode's inertia 
(e.g. lAerts et al.l 2010), normalised by the inertia of a ra- 
dial mode, interpolated to the same frequency. The kernels 
K™ l Ti (r) and Kfif p (r) are known functions which can be 
deduced from the (n, £) eigen mode of the reference mode l, 
via the variational principle (jGoueh fe Thompsonl Il991). 
Although the above relation can also be applied to other 
structural pairs, we will work with the pair (p, T±) since p is 
an obvious choice for finding the mean density, and Ti is ex- 
pected to vary little between different models. Throughout 
this paper, we will use the following definition for the rela- 
tive difference of a quantity /, a lthough we note tha t other 
definitions have been used (e.g. lAntia fe BasdfT994T ): 



Sf _ Jobs — /rcf 
/ /ref 



(2) 



With this definition, we have the relation / Q b s = /rcf(l + 
Sf/f). 

The mass difference, SM, between the reference model 
and the star is given by: 



SM = / 4nr z Spdr. 
'o 



(3) 



From this equation it is possible to deduce the relative dif- 
ference in mean density, p — 3Af/47ri? 3 : 



Sp 
P 



Attx 2 — — dx, 
PR P 



(4) 



where the quantity pn = M/R 3 is used for non- 
dimensionalisation. It is important to note that if Sp/p = 1 
then Sp/ p = 1, in other words: 



1 = 



Airx 2 — dx. 
PR 



(5) 



If we now assume that the radius of the reference model 
is different from that of the star, then the model can be 
rescaled through homology, i.e. a transformation which pre- 
serves all dimensionless quantities, so as to have the same 
radius as the star. Furthermore, the mean density needs 
to be kept the same, so as to avoid modifying the pulsa- 
tion frequencies when rescaling the model. Equation ([3]) 
will then give the relative mass difference between the star 
and the rescaled reference model. However, given their non- 
dimensional form and the fact that the density is preserved 
between the original and scaled models, Eqs. ([I} and (j3| re- 
main unchanged and can be applied directly to the original 
(unsealed) reference model, as long as x = 1 corresponds to 
both the radius of the star, and the radius of the reference 
model. The procedures which will be described below are 
based on these two equations, and therefore do not require 
explicitly scaling the reference model to the same radius 
as the star. They will consequently yield the mean density 
rather than the mass, although the later can be deduced if 
the radius is known. 



2.2. Linear inversions 

The procedures described in this paper then consist in 
searching for a linear combination of the relative frequency 
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differences that best reproduces the relative mean density 
variation, Sp/p: 



Spi, 



5ui 



(6) 



where the index i represents the pair (n,i). Using Eq. (|T|) 
to replace bvijvi then yields the following expression: 



= / K livg (x) — dx+ [ K cmss (x)^-dx 
Jo P Jo 1 1 



P 



where: 

K avg (x) = averaging kernel = ^ Cjif l p Tl (x), 



s (x) = cross-term kernel = \J CiKj< uf) (x). 



(J) 



(8) 
(9) 



In order for 5pi nv /p to be a good estimate of the relative 
mean density variation, the averaging kernel, A avg , needs 
to be as close as possible to 47rir 2 p/pR, as can be seen by 
comparing Eqs. (J7J and (@|. Furthermore, K CIOSS , which rep- 
resents the cross-talk coming from relative variations on the 
Ti profile, needs to be reduced as does the last term, which 
represents the surface effects. Hence, these kernels provide 
a natural way of evaluating an i nversion procedure (e.g. 
Christcnscn-Dalsg aard et al.lll990t ) . 



2.3. Various error bars 

The starting point for obtaining the error on the mean 
density estimate is to calculate the difference (5pi nv — 

£>Pexact)//0: 



Spinv - Spe 



K aV c — 47T— x 2 ) —dx 

(i \ PR J P 



^-cross t-, -f- ^ ^ 



£< 



(10) 



where the are the actual errors on the relative frequency 
differences resulting from measurement errors. If the true 
density and I\ profiles are known beforehand (for instance 
when testing out the methods on stellar models rather than 
true stars) then it is possible to quantify separately the dif- 
ferent contributions to the error using Eq. (fTO)) . Of course, 
in a realistic case, such profiles are not accessible. 

Using Eq. (|10j) . one can seek an upper bound on the 
error through the Cauchy-Schwartz inequality: 



Spinv — <5/?exact 



< 



A ave - Att-^x 2 



A", 



cross 1 1 2 



PR 

sr 1 





dp 


2 


P 



hi 



\ 




\ 




-(id 



where II/Hj = \/ '/„' P(x)dx. The difficulty, of course, is 
finding reasonable estimates for ||5p//9|| 2 , ||<5ri/ri|| 2 , and 
X^i [-fsurf (yi)IQi\- In their solar structural inversions, 



iRabello-Soares et al.l (|l999ft used the quantity ||A" avg — T|| 2 
as a measure of how well the averaging kernels matched 
specified target functions, T (which in their case were lo- 
calised around specific grid points), and ||-Kc ross ||2 to char- 
acterise how well cross-talk was being suppressed. 

It is also useful to look at the statistical error on the 
estimated density variation, resulting from frequency mea- 
surement error. Let us denote as <7i the la error bar on 
the relative frequency differences, Svijvi. The resulting la 
error bar on the estimated mean density variation is then 
given by: 



a Sp/p 



IE' 



(12) 



where we have assumed that the errors on the indivi d- 
ual frequencies are independent. The quantity \J^2 i c'j is 
known as the error magnification and is the amount by 
which the error is amplified if the o~i are uniform. 

It is important to realise that there are further sources 
of error besides those given in Eq. pop . These include non- 
linear effects where Eq. (fTJ) is no longer valid, non-adiabatic 
effects which are not included in the pulsation calculations 
nor in the variational principle (used to calculate the ker- 
nels), inaccuracies in the structural kernels which come 
from neglecting surface terms in the integration by parts 
(although these may be absorbed up to some extent in 
Si CiFsaxfiyi) / Qi) an d numerical inaccuracies. 

2.4. Homologous transformations 

Before going on to describe various inversion procedures, it 
is useful to point out simple properties which can be de- 
duced from homologous transformations. Let us suppose 
that the observed star can be deduced from the model 
through a homologous transformation, such that the rel- 
ative density variation is Sp/ p = e, where e is a small quan- 
tity. As is well known from dimensional analysis, the fre- 
quencies scale with y/ GM/R 3 . Therefore, to first order the 
relative frequency variation will be e/2. As a consequence, 
a necessary and sufficient condition to obtain the correct 
inversion result in the case of a homologous transformation 
is: 



£« 



(13) 



where we have made use of Eq. ©. In what follows, inver- 
sion procedures which satisfy this condition will be called 
"unbiased" . 

Furthermore, one can also deduce a simple property 
about the K P: y 1 kernels. Remembering that <5Fi/Li = 
in a homologous transformation and making use of Eq. ((TJ 
yields the following result: 



\* f K l pXi {x)dx. 



(14) 
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This relation is approximate because in deriving the ker- 
nel Kp^u there are several i ntegrations by part in whic h 
surface terms are neglected ([Gough fc Thompson! Il991h . 
Nonetheless, Eq. (fT4"]l is satisfied to a good degree of ac- 
curacy. 

Using Eqs. (fT4")l . and (O, it is possible to re-express 
Eq. (H31) as follows: 



K avg (x)dx : 
o 



Attx 2 —dx, 
PR 



(15) 



which is another way of expressing that an inversion pro- 
cedure is unbiased. 



3. Inversion procedures 

We will now describe three different methods for estimating 
the mean density of a star. The first method is an inver- 
sion procedure which yields directly inversion coefficients, 
Ci . The next two methods are non- linear scaling laws rather 
than inversion procedures. In order to put them in the same 
form as an inversion procedure, we linearise these laws and 
obtain associated coefficients c^. This then allows a system- 
atic comparison of the three methods using the averaging 
and cross-term kernels. Hence, although the second and 
third method are not initially inversion procedures, they 
will be treated and described as such in what follows. 



3.1. SOLA inversion 

The S OLA procedure (subtractive optim ally localised av- 
erages, iPiipers fc Thompson! [l992l 119941) consists in opti- 
mising K avg by minimising it's difference with a specified 
target function, T, and reducing the effects of the remain- 
ing terms. This is achieved by minimising the following cost 
function: 



J(Ci) 



{K avg (x) - T(x)} 2 dx + (3 I {K cmss (x)} 2 dx 



10 

A4„ rf 



J2 a ™Yl 

rn—l i 

tan 9 J^i (ci&iY 



Q, 



(16) 



where (er 2 ) = -k Y^iLi °f > N being the number of observed 
frequencies. In order to obtain the relative mean density 
variation, the appropriate target function is: 



T(x) = Attx 



2_P_ 

PR 



(17) 



The first term on the right hand side of Eq. ([T6]) is used to 
reduce the differences between K avg and T, and the next 
two terms minimise the effects of cross-talk from Ti and 
pollution from unknown near-surface effects. The functions 
ip m are a basis of slowly varying functions, such as the first 
Af sur f Legendre polynomials, and the CLyyi fire Lagrange mul- 
tipliers which are used to suppress these terms. The fourth 
term ensures that the inversion is not biased, A being a 
Lagrange multiplier. Finally, the last term is a regularisa- 
tion term which reduces crgp/p (see Eq. H"2"j) . The free pa- 
rameters in this approach are f3 which regulates between 
optimising K avg and minimising -Kcross, 8 which is used to 



adjust the amount of rcgularisation, and M sur f which con- 
trols the number of terms used in suppressing unknown 
surface effects. 

Using Eq. (fT7|) to define the target function avoids in- 
verting for the entire density variation profile, dp, and inte- 
grating it to obtain the mean density variation, which sim- 
plifies the calculations and reduces the number of free pa- 
rameters (such as the choice of a grid on which to carry out 
the inversion, or the width of the target functions at each 
grid point). Furthermore, the averaging and cross-term ker- 
nels are calculated for the inverted mean density variation, 
rather than for the density at specific grid points, which is 
more useful in ev aluating the in version errors. Finally, it is 
worth noting that lPiipersI (|1998l ) had already used a similar 
global approach to determine the sun's moment of inertia, 
and had found similar benefits. 



3.2. Large frequency separation 

An alternate approach is to use the well-known scaling law 
between the mean den sity and the large frequency separa- 
tion (e. g. |ulric3[l986|): 



(Av) oc yfp. 



(18) 



where (Av) is an average of the large frequency separa- 
tion, v n +i,l ~ v n.i- This law has been applied by many au- 
thors to obtain the mean density of a variety o f stars (e.g. 



Kieldsen fc Bedding! 119951 iMosser et all I2010D . Recently, 
White et al.l (|2011|) has investigated its accuracy for a grid 
of models and shown that the error is around 2 to 3.5 % for 
stars with T e g below 6700 K, relevant to solar-like oscilla- 
tions, but can reach 10% in the worst cases, when calibrated 
to solar values. 

Equation (TT81 stems from the fact that the frequencies, 
and consequently any linear combination thereof, also scale 
with ^fp for homologous transformations, as was pointed 
out in Sect. 12.41 This still approximately holds for non- 
homologous transformations, although some dependence on 
the precise stellar structure will appear. Asymptotic analy- 
sis shows that for high order p-modes, (Av) — 2 & 

i.e. the inverse of the time it takes for a sound wave to 
travel from the surface of a star to its centre and back 
again (|Vandakurovl 1 19671) . Given that sound waves travel 
more slowly near the surface, the large frequency separa- 
tion is sensitive to surface conditions. 

As such, this law is not directly comparable to the above 
procedure. However, one can apply this law in differential 
form, i.e. chose a reference model which is close to the star, 
calibrate the law using this model, and evaluate the mean 
density of the star. Given that the reference model is close 
to the star, one can simply perturb Eq. (fT8|) to relate the 
relative variation of the large frequency separation to the 
relative mean density variation: 



5(Av) 



5p 



(19) 



If we assume that the average large frequency separation 
can be expressed as a linear combination of the frequencies, 
(Av) = J2i d%Vi, the left hand side of Eq. (|T9)) becomes: 



S(Av) _ 2Y^di&Ui _ 2EidiVi S -% 



(Av) 



J2i diVi 



fc 7/. 



4 



D. R. Reese et al.: Estimating stellar mean density through seismic inversions 



where 



(21) 



j "3 "J 



Equation (|19p then takes on the form given in Eq. ([5]). 
Consequently, it is also possible to construct averaging and 
cross-term kernels for this approach, which can then be 
compared with those coming from the SOLA approach. By 
construction, this approach is unbiased and will therefore 
provide the correct answer for homologous transformations. 
One can also easily verify that J^i c i = 2. 

It important to realise that although the large frequency 
separation is typically applied to p-modes in the asymptotic 
regime, Eqs. (fT5)| and JT2J) have a more general application. 
Indeed, what is important is the y/p scaling which applies 
to all modes, not the fact that the large separation asymp- 
totically gives the inverse sound travel time. Hence, unless 
specified otherwise, we will typically use all of the available 
modes when applying Eqs. (fT5|) and ([TO)) . 

In what follows, we will use the following estimate for 
the large frequency separation in order to obtain the coef- 
ficients dj\ 



{Av 
where 
(Au(£)) 

W)) 



Egi fc(i)-M*))]Mi)-(rc(i))] 



Et e i [ni(<?)-(n(e))Y 



t=l 

N, 



N e 



(22) 

(23) 
(24) 
(25) 



Here, Ng is the number of modes with a given harmonic de- 
gree, £, and the rii(£) and vt{£) are their associated radial 
orders and frequencies, respectively. The above definition is 
a weighted average of the least-sq uare estimates of the m ean 
large frequency separation from iKieldsen et al~l (|2008l) for 
each I value. Using a least-squares approach has the advan- 
tage of not requiring consecutive radial orders to estimate 
the large frequency separation. 



3.3. The KBCD approach 



Recently, IKieldsen et al.l (|2008t ) proposed a method for cor- 
recting unknown surface effects which they treated as a 
power law. This method uses a combination of the large 
frequency separation and an average of the frequencies to 
estimate either the exponent of the power law or the mean 
density of the star. Havi ng established the valu e of the ex- 
ponent in the solar case, IKieldsen et al.l (|2008f) then go on 
to use this value to estimate the mean density of other 
stars. In what follows, we will refer to this method as the 
"KBCD approach" . Like the scaling law based on the large 
frequency separation, this approach is non-linear and there- 
fore not directly comparable to linear inversion methods. 
However, if one assumes that the reference model is close 
to the observed star, then it is possible to linearise this ap- 
proach as was done above. To do so, we start with Eq. (6) 



of IKieldsen et al.l (|2008lL subtract 1 from both sides, 
retain only first order terms: 



and 



ISp 
2j 



("> 



(Au) 



(26) 



where b is the exponent involved in the power law describ- 
ing the unknown surface effects, and (v) an average of the 
frequencies. Interesting limiting cases of the above equation 



• 6 = 0: this corresponds to fitting the large frequency 
separation and brings us back to the previous section, 

• b — > oo: this corresponds to fitting the average fre- 
quency, (y), directly. 

One will also note that b = 1 is a singular case which cor- 
respon ds to a degeneracy in the equations. IKieldsen et al.l 
(2008]) found b = 4.9 for the sun and consequently used this 
value for other stars. We will do likewise in what follows. 

In practise, we will use Eq. (1221) to obtain the large 
frequency separation, and a similar equation to obtain the 
average frequency: 



E| (Nt-i)W)) 

Y.iNt-1 



(27) 



thereby extending this method to non-radial modes. Once 
more, it is straightforward to see that this approach is un- 
biased. 



4. Non-linear extension 

The above procedures are linear and are therefore only 
adapted to reference models which are sufficiently close to 
the observed star. However, one may try to extend their 
applicability by pre-scaling the reference model so as to 
approach the mean density of the star and hopefully bring 
the relative differences into the linear regime. We therefore 
introduce a scale factor s, such that the mean density of 
the scaled reference model is s 2 p re f and the corresponding 
frequencies sv\ ct . With this rescaled model, the frequency 
shifts then take on the following expression: 



,rcf 



=- — +1 "I 



The inverted mean density becomes: 



Pinv(s) = /OrofS 2 < 1 + ^ ■ 



1 fSui 

S \ V4 



(28) 



(29) 



If one assumes that the linear inversion procedure is un- 
biased, then Eq. (f2Uf can be simplified as follows, using 
Eq. 



(30) 



Pinv(s) = Pre£ { -S 2 + S 




Equation (|30p is a 2 nd order equation with the following 
maximum value: 



Pmax — Pinv (^max) — Pref^ 



2 

max 5 



(31) 
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where 



1^ ^ 



(32) 



It is this maximum value which we will use as the best 
mean density estimate. Indeed, if one were to start with a 
given reference model, obtain a first mean density estimate 
through a given linear inversion, scale the reference model 
using this mean density, re-apply the inversion to obtain 
a second mean density, re-scale the reference model, and 
iterate till convergence, then the final mean density would 
be Pmax- Another way of looking at this is that s max is 
the scale factor for which the inversion yields no further 
correction to the mean density: 



1 fSvi 




- 1 













0. 



(33) 



It is then interesting to compare the result from a non- 
linear procedure to the p~max value obtained from the lin- 
earised version of the procedure. If we start with the scaling 
law between the large frequency separation and the mean 
density, then the p max value from the linearised procedure 
is: 



Prof 



1 



1 28 (Az/> 

2 (Ai/) 



Pref 



(Al/) 



(Ai 



' rcf 



(34) 



where we have used Eq. ©. This is identical to what is 
given directly by the original, non-linearised, scaling law, 
hence providing further justification for choosing p max as 
the mean density estimate. If we calculate p max for the lin- 
earised version of the KBCD approach, the result is: 



Mo 



Prof 



M rcf (A^) rcf 

6-1 



(35) 



do not depend on the scale factor. As such, they remain a 
useful diagnostic of the inversion procedure. However, when 
calculating the error as based on Eq. (fl"U|) . the relative dif- 
ference on the density profile needs to take into account the 
scaling by s max : 



scaled 



calcd 



(37) 



5. The solar case 



In order to illustrate the above methods, we apply them 
to the sun. A s et of 104 GOLF frequencies and associ- 
ated error bars (|Lazrek et al J 11997) were used as observa- 
tions. We only used the frequencies with a signal to noise 
ratio greater than 1, i.e. those which were not in brack- 
ets in Table I of lLazrek et all (|1997ft . The n and I val- 
ues of these modes are given in Table [TJ and the average 
of the error bars on these frequencie s is 0.694 /iHz. Model 
S ()Christensen-Dalsgaard et al1ll996[ ) was used as the ref- 
erence model. It is well known that there are surface ef- 
fects in the sun which are not accounted for in Model S. 
This causes the high order modes to differ in frequency 
from the observed values, as is illustrated in Fig. [TJ For 
the sake of comparison, the large frequency separation is 
also indicated on the fi gure. The pulsation modes w ere cal- 
culated with ADIPLS (IChristensen-Dalsgaardll2008l) . using 
an isothermal outer boundary condition, and the inversions 
were carried out using an improved (and corrected) version 
of InversionKitQ. 

Tab le 1. The n and I values of the 104 modes from 
lLazrek et ail (|1997t ). 

£ 1 2 3 4 5 = 

n 10-34 9-34 10-31 13-30 14-24 18, 19 



This time the result is different from what is given by 
the original non-linear procedure. However, this difference 
turns out to be minute. Indeed, Eq. ((35)) corresponds to 
the result which would be obtained had the unknown sur- 
face effects been represented as a power law in terms of the 
reference frequencies rather than the observed ones. Using 
this formulation, the exponent obtained in the solar case 
becomes 4.86 instead of 4.89, and mean density estimates 
are hardly affected (for b fixed at 4.9). Hence, in what fol- 
lows, we will use this formulation as representative of the 
KBCD approach. 

Finally, the associated la error bar on the estimated 
mean density is now given by: 



p IC {(T s 2 



.2^2 



(36) 



where we have assumed that ai -C 1. This is different from 
what is given in Eq. (|12[) due to the extra term s max . 

By construction, p max does not depend on the mean 
density of the initial reference model as opposed to a 
linear inversion result (this can also be shown explicitly 
through Eos. 1311 and l3"2"l). The same also applies to o"p maj[ . 
Furthermore, the various individual kernels, the averaging 
and cross-term kernels, and the target function remain un- 
changed since they are non-dimensional and consequently 



Before giving the inversion results, it is useful to dis- 
cuss what is the latest and hopefully most accurate value 
of the solar mean density. The mass and r adius used i n 
Model S correspond to the values given in 
Since t hen, the solar radius underwent a d ownward re- 
vision (jBrown fc Christensen-D alsaaard Il998f ) following a 
discrepancy with t he radius in ferred from solar f-modes 
(ISchou et aTJ Il997l lAntial Il998l). Further revisions can be 
found in Ha berreiter et al ] (12008ft who also point out the 
need for a generally agreed upon definition. The solar mass 
is determined by calculating the ratio GMq /G, where GM Q 
is known from planetary motions and G, the gravitational 
constant, from various experiments. The product GMq is 
known to a high degree of accuracy so the main source of 
uncertainty comes from the value of G. A summary of the 
latest values as well as those used in Model S can be found 
in Table [5] As can be seen in Table [5] (fifth row) , the latest 
mean density represents a 0.14 ± 0.09% increase over the 
value used in Model S. However, given that the pulsation 



1 This software can be downloaded form the Global 
Helioseismology HELAS website: 
http : / /helas . group . shef . ac . uk/ science/ 
inversions/InversionKit/ index . php 
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1.5 2.0 2.5 3.0 3.5 4.0 4.5 
v ohs (mHz) 

Figure 1. Differences between observed solar cyclic fre- 
quencies and those obtain ed from Mod e l S, u sing the I = 
— 5 frequency set from lLazrek et al.l (|1997f ) . The verti- 
cal grey lines indicate the la observational error bars. The 
large frequency separation is also indicated. 



pected from the latest solar parameters. This is likely to be 
caused by poorly modelled surface effects. However, we do 
note that the (Av) scaling is much more affected, which is 
not surprising given that the sound travel time integral is 
dominated by the near-surface contribution. The next col- 
umn gives the la error bar around Sp/p resulting from the 
measurement uncertainties on Sui/ui. Finally the last two 
columns give ||A" avg — T\\2 and 1 1 AT C ross 1 1 2 , which intervene 
in the upper bound on the error, as given in Eq. (fTTj) . 

As for most inversion methods, the results are a trade- 
off between minimising the effects of measurement errors 
and cross-talk from STi/Ti and optimising the averaging 
kernel. The parameters j3 — 10~ 8 , 8 — 10 -4 seem to be a 
good compromise. If /3 is increased to 10 _1 then the inver- 
sion wastes much effort trying to reduce the cross-talk from 
STi/Ti, which is expected to be small. Although || Across || 2 
is somewhat reduced, the inversion result is worse and asp/p 
has noticeably increased. If instead, 9 is reduced to 10~ 6 , 
the inversion result is slightly better as is || K avg — T|| 2 , but 
both asp/p and ||ifcroBs||2 take on large values. 



frequencies scale with y^GM/R 3 , it is not possible to dis- 
tinguish between a variation on the solar mass and a mod- 
ification of the gravitational constant from frequency data 
alone. Hence, it is more appropriate to compare Gp~Q be- 
tween Model S and the sun. As can be seen in the last row, 
the latest value represents a 0.16 ± 0.08% increase over the 
value from Model S. This corresponds to what one would 
expect from inversions if these were exact. Given the excel- 
lent agreement between GMq from Model S and the latest 
value (Table [21 first row), the main contribution to this 
difference comes from the discrepancy on the radius. 

Table 2. Solar mass and radius. 





Model S 


Latest values 


GM Q (lO^cm^s -2 ) 


1.32712445" 


1.32712440" 


G (10 -8 cni 3 g -1 s -2 J 


6.67232 


6.67384(80) c 


M Q (I0 33 g) 


1.989 


1.98855(24) d 


Rq (I0 10 cm) 


6.9599 


6.95613(171) e 


Pe (gem -3 ) 


1.4084 


1.4104(12) 


Gp Q (10- 8 s- 2 ) 


9.3975 


9.4128(70) 



Based on the product GMq using the values from the next 
two rows. ^ Based on ICoxl (|2000T ) and references therein. 

(c) 2010 CODATA recommended value, taken from 
http : //physics .nist . gov/cuu/Constants/ index.html 
^ Based on the ratio GMq /G using the values from the two 
previous rows. ^ Based on an a verage of the two values and 
quoted error bar in Table 3 of lHaberreiter et all (|2008l 1. 



Tabic [3] lists the results for SOLA inversions with var- 
ious parameter settings, and for a (Av) scaling. The cor- 
responding averaging and cross-term kernels are shown in 
Fig. [JJ The first two columns give the values of f3 and 9, 
when relevant. No surface corrections are included in the 
SOLA approach (i.e., M sur £ = in Eq. [T5]). The third col- 
umn gives the relative mean density variation (taking into 
account the non-linear generalisation described in Sect. [4]). 
As can be seen, all of the inversions predict a slight decrease 
of the mean density, as opposed to the 0.16% increase ex- 



5.1. Near-surface effects 

It is then interesting to look at different ways of correcting 
surface effects. Traditionally, these effects have been mod- 
elled as an additive term made up of a unknown function 
of frequency only, F sur f , normalised by Qi, the mode iner- 
tia divided by the i nertia of a radial mode interpolated to 
the s ame frequency (|Christense n-Dalsg aard fc Berthomieul 
11991ft . One way of dealing with this term in a SOLA inver- 
sion is to treat the function F sur f as a slowly varying func- 
tion of v, such as a linear combination of the first M sur f 
Legendre polynomials, and then t o constrain the inve rsion 
to cancel out any such function (jDappen et al.|[l99lT also 
see Eq. (TTBIO . The first row of Table 2] gives the inversion 
result and characteristics for Af sur f = 1, which amounts to 
treating F sur f as a constant. The corresponding averaging 
and cross-term kernels are given in the top panels of Fig. [31 
As can be seen, the results are quite poor. A likely reason 
for this behaviour is that the function F sur f is not well de- 
scribed by a constant. Instead, judging from Fig. [TJ and 
remembering that Qi is normalised in such a way as to be 
close to 1 for low £ values, F sur f starts at small values and 
increases with frequency. 

One way of trying to deal with this problem is to nor- 
malise the function F sur f in a different way, before treating 
it as a linear combination of slowly varying functions. Two 
alternatives are: 



Gsurf iyi) 

E,. ' 



and 



.-6+1/ 



(38) 



where we have used the notation G sur f and H sur { to dis- 
tinguish these functions from the original function F sm f, 
and where Ei is the mode inertia and b the exponent 
from the KBCD approach. The first possibility is based 
on the following reasoning. In the original description, the 
quantity Qi was introduced to replace Ei as it doesn't 
vary by seve ral orders magnitude over the p-mode fre - 
quency range ()Christensen-Dalsgaard fc Berthomieulll991ft . 
However, factoring out this strong variation from Ei also 
means introducing it into F SUT {. Therefore, normalising by 
Ei instead of Qi avoids this problem. The second possibil- 
ity is based on the idea that F snl {(vi) = Qi(6vi/h>i) SUI [ ~ 
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Table 3. Inversion results for the sun. 



p e 


Sp/p 




°~Sp/p 




l|AA' avg || 2 


\\K cross ||2 


10 -a 1Q -4 


-1.9 x 10" 


-3 


5.3 x 10 


-4 


0.32 


2.54 


10" 1 io~ 4 


-5.1 x 10" 


-3 


1.9 x 10 


-3 


0.35 


1.80 


10 -8 1Q -6 


-1.2 x 10" 


-3 


1.2 x 10 


-2 


0.31 


40.6 


(Au) scaling 


-1.2 x 10" 


-2 


4.1 x 10 


-4 


1.36 


2.77 



Notes. The first three rows correspond to SOLA inversions with different parameter settings and no surface corrections (M sur f = 0), 
whereas the last row corresponds to a { Av) scaling. The associated kernels are shown in Fig. O 




0.4 
Position 



1.0 



Figure 2. Averaging and cross-term kernels from solar inversions for various parameter settings. No surface corrections 
are included in the inversion (M sur f = 0). 



Table 4. Solar inversion results using different schemes for correcting surface effects. 



Description 


M surf 


8p/p 




°~Sp/p 




l|A^avg|| 2 


| -Across | 2 


Qnt normalisation 


1 


-5.6 x 10" 


■2 


2.2 x 10" 


•2 


1.84 


47.0 


E n £ normalisation 


1 


-1.0 x 10" 


-3 


1.2 x 10" 


-3 


0.50 


2.55 


v~ b+1 Q n t normalisation 


1 


-8.5 x 10" 


4 


6.2 x 10" 


-4 


0.40 


2.25 


v~ b+1 Q n i normalisation 


7 


-5.8 x 10" 


-2 


2.3 x 10" 


2 


1.86 


45.9 


Kieldsen et al. (2008) 




-1.9 x 10" 


-3 


5.5 x 10" 


5 


0.41 


2.03 



Notes. The parameters ft = 10 
in Fig. O 



6 = 10 4 were used for the SOLA inversions (first four rows). The associated kernels are shown 



.6-1 



av\ based on iKieldsen et al.l (|2008l ). Hence by dividing 
both sides by v"^ 1 , one is left with a function H sul f which 
is slowly varying. 

Rows 2-4 of Table 0] give the results from applying these 
two alternative normalisations, and the corresponding ker- 
nels are shown in Fig. [3] In the first two cases, only one (con- 
stant) term is used to correct surface effects. The results 



are slightly better than when no surface corrections are in- 
cluded, with a slight preference for the v~ h+1 Qi normalisa- 
tion. At the same time, both the averaging kernel and &5p/p 
are slightly worse. In last case, M sur f = 7 in accordance with 
the op timal number of terms found in lRabello-Soares et al.l 
(|1999D . However, this leads to poor results. The main dif- 
ference between the inversions carried out here and those 
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carried out in Rabcllo- Soares et al.l (|l999ft is the number of 
modes involved in the inversion. Given the limited number 
used here (which still remains generous in comparison with 
stars other than the sun), one cannot afford to over-correct 
for surface effects. However, even when no surface correc- 
tions are explicitly included, Table [3] shows that a SOLA 
inversion is less affected by surface effects than a (Az/) scal- 
ing, probably because as K avg approaches the target func- 
tion T, its surface amplitude is reduced. 

Another way to try to limit surface effects is to perform 
a SOLA inversion using o nly kernels based on frequency 
separation ratios. Indeed, iRoxburgh fc Vorontsovl $2003) 
showed that these ratios are insensitive to surface condi- 
tions through analytical developments based on internal 
phase shifts and through numerical simulations. Likewise, 
(|Oti Floranes et al1l2005l ) showed that the associated ker- 
nels have a small amplitude near the surface. However, they 
also pointed out that both the numerator and denomina- 
tor of a frequency separation ratio scale with -J GM/R 3 , 
meaning that the ratio is insensitive to the mean density of 
the underlying model. As a result, the Cj coefficients associ- 
ated with an inversion based on these kernels add up to 0, 
making it impossible to do an unbiased inversion. Likewise, 
it can be shown using Eq. (fT4")l that the integrals of the 
kernels K la ^° is approximately 0. 

Finally, the last row of Table @] gives the results from 
applying the KBCD approach. The associated averaging 
and cross-term kernels are displayed in the lower panels 
of Fig. |31 The mean density estimate and crgp/p compare 
favourably with the other inversion results and the asso- 
ciated kernels are well behaved, thus explaining why this 
method works. This approach is found to be superior to 
the simple scaling law based on the large frequency separa- 
tion, as can be seen both in terms of the results it produces 
and in terms of the kernels. The likely reason for this is that 
the averaged frequency, which intervenes in this approach, 
uses a more uniform weighting on the individual frequencies 
whereas the large frequency separation favours frequencies 
at either end of the frequency range, thus amplifying surface 
effects which are strongest at high radial orders. Overall, 
the best results are obtained with the v~ h+1 Qi normali- 
sation, using one term for surface effects, and the KBCD 
approach (3 rd and 5 th line of Table 0}. 



with some evolutionary tracks, the dashed parts corre- 
sponding to the prc-main-scqucncc. Three additional mod- 
els, nicknamed "Models A, A' and B", were selected near 
the centre of this box to play the role of observed stars. 
These models are also shown in Fig. [U and Table [5] gives 
some of the characteristics of Models A and B. The fre- 
quencies of these models are shown in an echelle diagram 
in Fig. [3] In all the cases tested here, the la error bar on 
the frequencies was set at 0.3, /HLz, although no errors were 
added to the frequencies. The goal of this study is then to 
see how well the density of these three models can be re- 
produced by applying the various inversion procedures to 
the models from the grid. 
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o 
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4.50 








4.45 




it- r,~ r, rr,7,\r,~ 7 





5100 5200 5300 5400 5500 
Teff (K) 



Figure 4. T c g - log(g) diagram showing the position of 
Model A (0), Model B (□), the reference models from grid 
B (+), and various pre-main-sequence and main-sequence 
tracks, represented by dashed and continuous line respec- 
tively. 



6. Test cases with a grid of models 

In order to do a more systematic study, we apply the dif- 
ferent inversion procedures to a grid of 93 main-sequence 
and pre-main-sequence models. These models were down- 
loaded from grid B of the CoRoT-ESTA-HELAS web- 
sitc0, and were chosen within the (T e g, log(g)) box (5291 ± 
100 K, 4.563 ± 0.1 dex), the error bars being representative 
of what can be expected from spectroscopic observations. 
Their masses range from 0.80 Mq to 0.92 Mq and their ages 
from 28 Myrs to 17.6 Gyrs, after their appearance on the 
birthline. A description of the physical ingredients, chem- 
ical compositi on and initial conditi ons used in these mod- 
els is given in lMarques et a l. (2008). In particular, the ini- 
tial chemical composition is (X,Y,Z) = (0.70,0.28,0.02), 
where the m etal abundances correspo nd to the solar mix- 
ture given in IGrevesse fc Noels! (11993ft . Figure H shows the 
positions of these models in a (T c s , log(g)) diagram, along 

2 http : //www. astro .up.pt/corot/models/cesam/ 



Table 5. Characteristics of Models A and B. 





Model A 


Model B 


9 (gem 3 ) 


2.2895 


2.3079 


Mass (Mq) 


0.900 


0.920 


Radius (Rq) 


0.821 


0.825 


Age (Gyr) 


1.492 


2.231 


log(s) (dex) 


4.563 


4.569 


Teff (K) 


5291 


5291 



Notes. Here, Mq and Rq take on the values from Model S. 



6.1. Model A 

Model A represents a 0.9 Mq main-sequence star and was 
downloaded from grid A of the CoRoT-ESTA-HELAS web- 
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0.4 0.6 
Position, r/R 



0.4 0.6 
Position, r/R 



Figure 3. 

shows the 



Averaging and cross-term kernels from solar inversions for different treatments of surface effects. The last row 
averaging and cross-term kernels from the KBCD approach rather than from a SOLA inversion. 



sitc0- Models from grid A are started off as fully convective 
spheres, instead of using a birthline as in grid B. Apart from 
this difference, this model has the same chemical composi- 
tion and physical ingredients as the models in the grid. 

Each one of the models from the grid was used as a 
reference model from which to apply the different inver- 
sion techniques to the "observed" frequencies from Model 
A. Figure |6] shows the results from these inversions for a set 
of 33 modes with radial orders n = 15 — 25 and harmonic 
degrees £ = — 2. The parameter j3 was set to 10~ 8 for all 

3 The specific URL to Model A is 

http : //www. astro .up.pt/corot/models/cesam/ 
A/data/O . 90/ms/star_1492 . 002-ad . osc . gz 



of the SOLA inversions. Although the mean densities of the 
reference models differ from that of Model A by up to 30 % 
(see the rc-axis in figures), the inversion results are within 
0.3 to 0.6% of the true mean density except for the worst 
cases where the departure can exceed 1.0%, which shows 
the sensitivity of these modes to density. A more detailed 
look at the figures shows that overall, the KBCD approach 
and SOLA inversions produce the best results, provided 
the parameter 9 is not too small. The two are very simi- 
lar if M sur f = 1, 8 = 10 -2 and a v~ h+1 Q n i normalisation 
is used for surface effects. The large frequency separation, 
for the most part, produces results which are worse. One 
may also notice that in most cases, the difference between 
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Figure 5. Echelle diagram with the frequencies from 
Models A, A' and B. The £ = 2 and £ = ridges are to 
the left, whereas the £ — 1 ridge is to the right. Note: the 
large frequency separation indicated in the lower right cor- 
ner is the value used to construct this hgure, not necessarily 
the exact value for the three models. 



the inverted and true mean densities is much larger than 
the la error bar from the inversions. This behaviour can 
be expected, since the la error bar from the inversion is 
only based on the frequency measurement errors and does 
not take into account the other sources of error mentioned 
in Sect. 12. 31 such as differences between the averaging ker- 
nel and the target function or cross-talk from the STi/Ti 
profile. 

As can be seen, reducing 8 mainly seems to increase the 
la error bar in SOLA inversions without improving the re- 
sults. This can be understood by looking at the trade-off 
curves illustrated in Fig. Indeed, reducing 8 only has a 
small effect on ||ifavg — T|| 2 while a^p/p increases substan- 
tially. Furthermore, the same mechanism which amplifies 
observational errors can also amplify errors which are not 
accounted for in Eq. (TIT)]) , such as non-linear effects not in- 
cluded in Eq. (fTJ) or numerical errors on kernel calculations. 
Based on Figs. [5] and [JJ the value 8 = 10~ 2 seems to be a 
good compromise between the two effects. 

Furthermore, given that the true structure of Model A 
is known, it is possible to analyse the different sources of 
error based on Eq. (jlOl) . These are displayed in Fig. [8] for 
SOLA inversions. The dominant contribution comes from 
a mismatch between the averaging kernel and the target 
function and is represented by a solid line. The dotted line 
represents the cross-talk from the STi/Ti profile and is the 
smallest contribution. Lastly, the dashed curve represents 



0.0020 
1 0.0015 



0.0010 
0.0005 
0.0000 
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0.46 0.41 



Higher 8 



3 0.50 0.52 

lK avg -Tll 2 



Figure 7. Superimposed trade-off curves for each of the 
93 models from the grid. The parameters are M sur f = 0, 
/3 = 1CT 8 and 8 = 1CT 6 , lfT 4 , 1CT 2 , 1, 1.5. The grey "x" 
marks the locations of intermediate 8 values. 



remaining effects such as surface contributions, measure- 
ment errors on the frequencies and effects not included in 
Eq. (flU)) . such as non- linear effects not accounted for in 
Eq. (TT]). Of these three, the non-linear effects are likely to 
be the strongest because the data is error free (except for 
numerical round-off errors) and no attempt was made to 
simulate surface effects in Model A. As can be seen, the 
non-linear effects are not negligible compared to the other 
sources of error. This is not entirely surprising given that 
the only criteria used for selecting the models from the grid 
was the (T e s, log(g)) error box, and consequently shows the 
limits of Eq. ((T|). It must also be noted that the errors in 
Fig. [8] tend to compensate each other. However, one cannot 
always expect such fortuitous cancellation to arise in other 
situations. 
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Figure 8. Different sources of error based on Eq. (jTOj) for 
SOLA inversions. See text for details (in particular for 
"Other"). The SOLA parameters are /3 = 10" 8 , 8 = 10~ 2 
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Figure 6. Inversion results for Model A, using the modes with n — 15 — 25 and £ = — 2. The x axis gives the mean 
density of the reference models divided by that of Model A, whereas the y axis is the mean densities obtained from the 
various inversion procedures, normalised by that of Model A. For each procedure, there are three lines corresponding 
to the result and the la error bar around it. The horizontal triple-dot-dashed line shows the correct mean density. A 
v~ b+1 Q n i normalisation has been used when surface corrections are included in SOLA inversions. 



Figure [9] shows what happens if the non-linear extension 
described in Sect. 0]is not applied to the inversion proce- 
dure. As expected, inversion results are substantially worse 
when the mean density of the reference is far from the true 
mean density. 

Finally, the wiggles observed in Fig. [6j especially for 
reference models with a lower mean densities, can be par- 
tially explained by the fact that the grid contains a mixture 
of pre-main-sequence and main-sequence models. Indeed, 
when the two types of models are separated as in Fig. ITUI 
two separate trends can be observed which likely reflect 
their differences in structure. 

6.2. Near-surface effects (Model A') 

Model A' is nearly identical to Model A, except for an ad- 
hoc 50 % increase of the density in the outer 1 % of the 
model, the transition being represented by a hyperbolic tan- 
gent function. The pressure was then recalculated using the 
hydrostatic equation, thus resulting in a very slight offset 
from the original profile. The difference in mean density be- 
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Figure 9. Same as Fig. [5] (right, middle panel), except that 
the non-linear extension described in Sect. 0] has not been 
applied to the inversion procedures. 



tween the two models is negligible, but the frequencies are 
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Figure 10. Same as Fig. [5] (left, middle panel), except that 
the pre- main-sequence models (left panel) and the main- 
sequence models (right panel) have been separated. 



noticeably different (see Fig. [S]) due to the fact that sound 
waves spend more time in the outer layers of a star. 

Figure [TT] shows the inversion results for this model, 
using the same set of modes as in the previous case. 
Overall, the results are worse than in the previous case. 
Furthermore, the scaling law based on the large frequency 
separation yields results which are substantially offset from 
the true mean density, when compared to the other pro- 
cedures. This could be expected given the fact that the 
large frequency separation, which asymptotically gives the 
travel time of a sound wave from centre to surface, is very 
sensitive to surface conditions. This high sensitivity to sur- 
face conditions is also reflected in the associated averaging 
kernels which tend to have large amplitudes near the sur- 
face. The SOLA inversions without surface corrections also 
seem to be affected, but to a lesser degree. In particular, 
it produces worse results than the KBCD approach for the 
higher density models, which is the opposite from the pre- 
vious case. If surface corrections are included in the SOLA 
approach, then these produce slightly better results than 
the KBCD approach. Once more, reducing 9 does not im- 
prove the SOLA inversions but only serves to increase the 
lcr error bar. 

6.3. Uncertainties in the stellar physics (Model B) 

Model B is different in a number of ways from Model 
A and also from the models in the grid. It has a differ- 
ent chemical composition and includes the effects of dif- 
fusion and rotational mixing, thereby providing a way to 
test the effects of unknown physics on mean density inver- 
sions. Specifically, the chemical abundances are (X, Y, Z) = 
(0.735,0.250, 0.015) and follow the more recent solar mix- 
ture given in lAsplund et all (|2005f ). 

The results are shown in Fig. Q21 using once more the 
same set of modes. These results are in fact similar to those 
obtained for Model A and shows that all of the inversion 
procedures presented here seem to be fairly robust to un- 
known physics. 

6.4. The effects of mode misidentification 

We then investigate the effects of misidentifying the radial 
orders of the modes. Figure Q2] (left panels) shows mean 
density estimates for Model A, using the same set of ob- 
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Figure 12. Same as Fig. [6] except for Model B. 

served modes as previously, but identified as n = 16 — 26 
(upper panel) and n — 14 — 24 (lower panel). As can be 
seen in these figures, the scaling law based on the large 
frequency separation provides the best results, whereas the 
other methods are off by 10%. A simple explanation is that 
the relative change over one radial order of Ais is smaller 
than that of v, given that Ai/ goes to a constant when n 
goes to infinity, whereas v is roughly proportional to n. 

Another way of misidentifying modes is by labelling the 
£ = modes as I — 1 and vice-versa {e.g. iBenomar et al.l 
2009a). This, however, also leads to detecting a different set 
of £ — 2 frequencies given that these will always be close 
to the £ = frequencies. In order to investigate this sce- 
nario, we swap the £ = and 1 identifications and create 
a fictitious set of £ = 2 frequencies in such a way as to 
preserve the original small frequency separation, averaged 
over adjacent radial orders. Two possible ways of assigning 
radial orders exist and are given in Table [SJ The correspond- 
ing mean density estimates are given in the right panels of 
Fig. Q2J Once more, the large frequency separation gives 
the best results, the other methods being off by about 5%, 
i.e. about half the amount as in the previous case when 
only n was misidentified. One way of interpreting this is to 
view this as a ±1/2 offset on the radial order, and then to 
apply the same argument as previously. 

Table 6. Different ways of assigning radial orders when 
misidentifying the £ value. 



True id. Upper n Lower n 

£ = 1, n = 15 - 25 £ = 0, n = 16 - 26 £ = 0, n = 15 - 25 

£ = 0, n = 15-25 l=\, n = 15-25 I = 1, n = 14 - 24 

Fictitious t = 2, n = 15 - 25 £ = 2, n = 14 - 24 



These observations provide a way of identifying modes: 
one can calculate the mean density using the large fre- 
quency separation and another method such as the KBCD 
approach or a SOLA inversion and search for the values 
of £ and n for which the two methods agree to within 1% 
or so. This procedure, although different from the methods 
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Figure 11. Same as Fig. [6] except for Model A'. 



descri bed in iBedding fc Kieldsenl (|2010ft and IWhite et al] 
(|201ltl , is based on the same principles . Indeed, the method 
presen ted here and the approach from IBedding fc Kieldsenl 
(|2010() both amount to scaling a model (or a comparison 
star) so that the large frequency separation agrees with 
that of the observed star, then adjusting the identification 
so that th e offset betw e en th e two sets of frequencies is 
minimised. IWhite et al.l (|201ll ) achieves the same result by 
directly comparing their model-based values for the large 
frequency separation and e, the additive offset from the 
asymptotic formula for p-mode frequencies, to those from 
observed stars. Of course, surface effects acting upon the 
large frequency separation may lead to an offset even for 
the correct identification, although we expect such an off- 
set to be smaller than what is shown in Fig. 111! which was 
calculated fo r a ra ther drastic modification of the surface. 
IWhite et al.l (|201ll ) have nonetheless pointed out situations 
where e can vary rapidly thereby making mode identifica- 
tion more difficult. 



6.5. Kernel mismatch 

It it also interesting to test the various inversion techniques 
on different sets of modes. We therefore consider the follow- 
ing 2 sets: 



• Set 1: n = 2 - 10 for I = 0, 2, 3, n = 5 - 10 for 1 = 1, 
n = 1 - 10 for £ = 4. 

• Set 2: The same as Set 1, without the n = 1, £ = 4 
mode. 



The mean density estimates are shown in Fig. 1141 As can 
be seen by comparing the two figures, including the n = 1, 
I — 4 mode has a drastic effect on the SOLA inversion re- 
sult for one of the models, which we shall call the "worst 
model". The explanation for this is quite simple. Figure IT51 
shows the kernels associated with this mode for both the 
worst model and for Model A. Clearly, the kernels in the 
worst model are dominated by a g-mode behaviour (even 
if the eigenfunctions still seem to be dominated by a p- 
mode character), whereas those of Model A are typical of 
p-modes. Such a difference on the kernels inevitably leads 
to a large discrepancy between the actual frequency shift for 
that particular mode, and the frequency shift which would 
be deduced by applying the linear approximation given in 
Eq. ([T]). This, of course, then affects any mean density esti- 
mates, particularly a SOLA inversion which explicitly op- 
timises the averaging and cross-term kernels using the in- 
dividual kernels. 
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Figure 13. Inversion results for Model A, using the same modes as in Fig. [5] but with various erroneous mode identifi- 
cations, as indicated in the panels. In all of the panels, j3 = 10~ 8 , 8 — 10~ 4 and M sur f = 0. 
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Figure 15. Kernels of n = 1, i = 4 mode for Model A (lower panels) and the worst model (upper panels). 



7. Observed stars 

7.1. a Cen B 

Having tested the inversion techniques on stellar models 
in a number of ways, we now try them out on observed 
stars. A particularly interesting example is the nearby Kl 
main-sequence star a Cen B. Its close proximity allows 
precise interferometr ic measurements of its angular radius 
(jKervella et al.1120031 ). When combined with a reprocessed 
Hipparchos parallax (Sodcrhiclm 1999), this yields a ra- 
dius R = 0.863 ± 0.005 R Q . Given that a Cen B is part 
of a multiple system, it is possible to determine it s mass 
through orbital parameters. Hence, iPourbaix et al.l (|2002l ) 



obtained a mass of M = 0.9 34 ± 0.0061 Mq using the par- 
allax from lSoderhielml (|1999T) . This then yields a mean den- 
sity of 2.046 ± 0.058 gcm~ 3 . A time series of radial veloci- 
ties was obtained bv iKieldsen et al.l (|2005h using the Very 
Large Telescope (VLT) and the Anglo- Australian Telescope 
(AAT). A Fourier analysis of this time series yielded a set 
of 37 frequencies with harmonic degrees £ = to 3, and 
radial orders, n, between 17 and 32. 

Given that the mass of a Cen B is close to the mass 
range covered by the grid from the previous section, we 
simply reused the models from this grid as reference models 
from which to estimate the mean density of a Cen B. The 
results are shown in Fig. 1161 These turn out to be similar to 
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(|2002h and the radius from iKervella et all (j2003h lead to a 
mean density of 2.046 ± 0.058 g cm -3 . The results from the 
KBCD and SOLA approaches are 2 to 3 % lower than this 
value, but remain within the error bars. They agree with the 
range of values, delimited by t he triple-dot-da s hed h orizon- 
tal lines in Fig. (TBI obtained in lKieldsen et al.l ()2008f ) where 
the KBCD approach is developed. The seismic study in 
lEggenberger et all (|2004[ ) leads to a mean density of 1.998 
gem -3 which also falls within this range, although they did 
increase the error bar on the radius to reach this value (the 
error bar on the mass seems to have been ignored in their 
study). On the other hand, the {Av) scaling leads to val- 
ues below all of these previous studies and which mostly 
lie outside the error bars of the non-seismic mean density. 
Hence, using the KBCD or SOLA approaches leads to bet- 
ter agreement with the non-seismic constraints and tends to 
confirm the idea that the (Av) scaling is affected by surface 
effects. 



Figure 14. Inversion results for Model A for two similar 
sets of modes (see text for details). The parameters used in 
these inversions are /? = 10~ 8 , 9 — 10~ 4 and M sur f = 0. 



the results obtained in the previous sections. The tendencies 
are the same, as is the spread in inverted mean densities, ex- 
cept for the results from the {Av) scaling where the spread 
is somewhat larger. Nonetheless, a noticeable difference ap- 
pears between the estimates from the {Av) scaling and the 
two other methods. This suggests that surface effects play 
an important role in this star. 



Table 7. Mean density of a Cen B. 



Source 




P 








(g cm 


- 3 ) 


Non-seismic™ 




1.988 - 


2.104 


Eggcnbcrgcr ct 


al. (2004) 


1.998 


Kieldsen et al. 


(2008) 


1.998 - 


2.011 


SOLA (M surf = 


0) b 


1.987 - 


2.006 


SOLA (M surf = 


1)" 


1.997 - 


2.012 


{Av) scaling 6 




1.952 - 


1.994 



Based on the parallax from Sodcrhielm] ll 19991) . the 
interferometric observations from IKervella et all (2003) and the 
orbital parameters from lPourbaix et al.l l|2002l ). This study. 
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Figure 16. Mean density estimates for a Cen B. The SOLA 
parameters are (3 = 10~ 8 and 8 = 10~ 2 , and Af sur f is in- 
dicated in the panels. The vertical triple-dot-dashed line 
corre sponds to the minimu m and maximum mean densities 
from lKieldsen et al.l (|2008t) . 



It is then interesting to see how these densities compare 
with previous studies. Table [7] gives a summary of previ- 
ous and current estimates of the mean d ensity of a Cen 
B. As was mentioned above, the mass from lPourbaix et al.l 



7.2. HD 49933 

We now turn our attention to stars observed by CoRoT. A 
particularly interesting example is the F5 main-sequence 
star HD 49933 which has been observed in two sepa- 
rate CoRoT runs, thereby yielding a rich frequency spec- 
trum. In what follows, we will use the set of 51 frequen- 
cies and associated err or bars and identifications given in 
(|Benomar et al.l [200 9b]). which takes into account the two 
observational runs. An alternative identification had been 
favoured before data from the secon d observational run was 
available (Appour chaux et al.ll200 8. e.g.). The grid of mod- 
els used in the previous sections is not suitable for this star. 
We therefore downloaded a set of 42 pre-main-sequence 
and main-sequence models from the CoRoT-ESTA-HELAS 
website, which lie in the error box T c tf — 6665 ± 215 K and 
lo g(I//L w ) = 0.54 ± 0.0 2 dex, based on the values quoted 
in lBenomar et all (|2010|) and references therein. It must be 
noted that no attention was paid to possible differences in 
chemical composition between these models and HD 49933. 

The mean density estimates are given in Fig. 1171 They 
fall in the interval 0.580 ±0.015 gem" 3 which corresponds 
to a ±2.6% range around the mid-value. This interval is 
on the lower end of the range of values covered by the 
reference (or input) mean densities. Furthermore, the es- 
timated mean density slightly increases with the reference 
mean density, in much the same way as in the previous 
examples except on a larger scale. If we limit ourselves to 
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Figure 17. Mean density estimates for HD 49933. The 
SOLA parameters are ft = 10 -8 and 9 = 10~ 2 . The hori- 
zontal triple-dot-dashed line s represent the solutio ns from 
ICreevev L Bazotl (|201ll) and IKallinger et all (|2010l ). in the 
same order as in Table [5] 

reference models with a closer mean density to the results, 
this would favour the lower values, but this reasoning must 
be taken with some caution since the inversions use the 
non-linear extension and are therefore insensitive to the 
mean density of the input models. Various other groups 
have performed detailed asteroseismic interpretation of his 
star and have come up with best fitting models using a va- 
riety of techniques. Here, we will only focus on studies for 
which both t he mass and radius of the best fitting model 
are provided. iKallinger et al.l <|2010T) searched for a best fit- 
ting model in a dense grid using a \ 2 criteria, whereas 
ICreevev fc Bazotl (|201lD found two solutions using both a 
Levenberg-Marquardt and a Markov Chain Monte Carlo 
approach. Their results are summarised in Table [5] and rep- 
resented as triple-dot-dashed horizontal lines in Fig.[T7J The 
first solution uses a mode identification in which the I val- 
ues are swapped. Its mean density is in the upper half of 
the range of values found here. The two other solutions 
use the same identification as the one used here. They have 
very similar mean densities in spite of having different radii 
and masses, and are at the lower end of the interval of so- 
lutions found here. Although further from our mid-value, 
these solutions are located where the mean densities are 
least modified by the inversions. 



Table 8. Parameters of best fitting models for HD49933 
from other studies. 



Source 


M 


R 


P 




(M e ) 


(Re) 


(gem" 3 ) 


Creevev & Bazot (2011). l a 


1.12 


1.39 


0.588 


Creevev & Bazot (2011). 2 


1.20 


1.44 


0.567 


Kallinger et al. (2010) 


1.325 


1.49 


0.565 



^ This solution is based on a swapped i identification. 



Given the difficulties in labelling the modes prior to 
the second CoRoT run, we tested out the swapped I iden- 



tification, using only the i — and 1 modes. The re- 
sults are shown in Fig. [18] and are very similar to the 
results presented in Fig. [T3] Clearly, the models used 
in this study favour the currently accepted identifica- 
tion. In par ticular, this agrees w ith the identification 
favo ured by Kallinger et al. (Hp iO) and with the results 
from lBedding fc Kieldsenl (|2010t) who compared HD 49933 
to other observed stars in which the identification is known. 
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Figure 18. Mean density estimates for HD 49933 us- 
ing swapped I identifications. The SOLA parameters are 
P = 1Q-8, = 1Q- 2 anc j M sur f = 1. The so l utions 
from ICreevev fc Bazotl (|201lD and IKallinger et"ai1 (|2010l) 
are represented as triple-dot-dashed horizontal lines. 



Finally, it is worth noting that surface effects have not 
manifested themselves in the results presented here: all 
three methods yield similar results for the correct identi- 
fication. This suggests that surface effects are less impor- 
tant in this star co mpared to other effect s deeper within the 
star. Interestingly, Bedding et al. (2010) recently reported 
on an asteroseismic investigation of Procyon, also of spec- 
tral type F5. They showed that the KBCD approach with a 
solar exponent, b = 4.9, did not provide a good description 
of the differences between the frequencies of various best 
fitting models and those of the star. Instead, these were 
better described as a constant offset, which corresponds to 
6 = 0. Given that these differences do not go to zero as the 
frequency is reduced, they also interpreted this as signify- 
ing that the structural differences between their best fitting 
models and Procyon extend beyond the surface layers. 

7.3. HD 49385 

The star HD 49385 is an evolved GO type star which 
was obser ved by CoRoT f or 137 days during its second 
long run. iDeheuvels et al.l (|2010l) were able to detect 27 
I = — 2 frequencies as well as what seems to be mixed 
modes and some unconfirmed I = 3 frequencies. In what 
follows, we will not use the I = 3 modes given their 
low signal to noise ratio, and will postpone dealing with 
the mixed modes to the end of this section, since these 
vary rapidly from one model to another, thereby limiting 
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the validity of Eq. ((T|). A new set of 18 models within 
the (r Bff ,loK(q)) = (6 095 ± 65K, 4.00 ± 0.06) error box 
(|Deheuvels et al.ll2010t ) was downloaded from the CoRoT- 
ESTA-HELAS website. Given the position of this box in 
the HR diagram, these models correspond to either pre or 
post-main-sequence starfl 

The mean density estimates are given in Fig. 1191 
The horizontal triple-dot-dashed lines are the mean den- 
sities corresponding to th e four families of solutions in 
iDeheuvels fc Michel! (|201lh . In the upper panel, the results 
are particularly poor. This failure is caused by differences 
in mode labelling between the pre and post-main-sequence 
models due to the presence of mixed modes in the lat- 
ter. The inversion results were based on an identification 
which corresponds t o some of the best-fitting models from 
IDeheuvels fc Michell (|201lf ) who applied Newton's method 
to a grid of models specifically designed to handle mixed 
modes. Such an identification is ill-adapted for pre-main- 
sequence models and to a lesser extent for some of the 
post-main-sequence models, due to the jumps it introduces 
in the n values. Consequently, the results from the five pre- 
main-sequence models, indicated by the circles, stand out 
as being the worst. In order to circumvent these difficulties, 
a second inversion was carried out, only using the I = 
modes. This time, the results a re far better and are within 
2.5% of the mean densities from lDeheuvels fe Michell ()2011l ) 
for the SOLA inversions and for the KBCD approach, and 
3.5% for the scaling law based on the large frequency sep- 
aration. Paradoxically, the pre- main- s equence models now 
give t he results which are closest to IDeheuvels fc Michell 
(|2011h . Once more, the scaling law is somewhat different 
than the two other approaches which is probably indicative 
of surface effects. 

In order to investigate the influence of mixed modes 
on SOLA mean d ensity inversions, we use a low over- 
shoot model f rom Deheuvels fe Michell (|201lh with the 
lAsplund et al.l (|2005l) mixture. Inversions were carried out 
for three different sets of modes: 



• Set 1 

• Set 2 

• Set 3 



26 p-modes with £ = — 2 

Set 1 + the (n =1A,£= 1) p-mode 

Set 1 + the (n = 10, 1 = 1) mixed mod<Q 



Two different sets of SOLA parameters were used: 

• Case 1: j3 = 10~ 8 , 6 = 10~ 2 and M surf = 1. These 
parameters were chosen in accordance with typical in- 
version parameters (the number of modes was sufficient 
to allow M sur f = 1). 

• Case 2: j3 = 0, 9 = 0, M sur f = 0. These parameters 
correspond to optimising K avg only. 

The results are summarised in Table H] As can be seen 
in case 1, adding a p-mode or a mixed mode has a similar 
impact on both the averaging kernel and the cross-term 
kernel. The reduction of the statistical error on the inverted 



4 As was argued in Deheuvels & Michel (2011), the presence 
of mixed modes shows that this is a post-main-sequence star. 
Hence, one would want to limit the inversions to this type of star. 
However, for the sake of comparison, we will also use pre-main- 
sequence stars and see what impact this has on the inversions. 

5 We note that there there is an offset of 1 between the ra- 
dial o rders used here and the ones used in Deheuvels & Michel 
(|201lT ) . This is merely a difference in mode labelling and conse- 
quently does not affect inversion results. 
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Figure 19. Mean density estimates for HD 49385. The 
SOLA parameters are /3 = 10~ 8 , 6 = 10 -2 and M surf = 0. 
The horizontal triple-do t-dashed lines corres p ond t o the 
best fitting models from IDeheuvels fe Michel (|201lh . The 
circled results (upper panel) were obtained for pre-main- 
sequence models (for clarity, these are not circled in the 
lower panel). 

mean density, ap, is larger for the mixed mode than for the 
p-mode, indicating that the mixed mode may be yielding 
more information on the mean density. This seems to be 
confirmed in case 2, where only the averaging kernel is being 
optimised. Indeed, the averaging kernel including the mixed 
mode is closer to the target function than the averaging 
kernel which uses the extra p-mode. However, the extra 
benefit from the mixed mode remains small and is probably 
of little use because the domain for which Eq. [T] is valid is 
likely to be much smaller for these modes than for ordinary 
p-modes. In some cases, non-linear effects on such modes 
can even cause a substantial error on the mean density as 
was shown in Sect. 16.51 

Table 9. SOLA inversion results for HD 49385. 



Mode set 


Pinv 

(gem -3 ) 


dp 

(gem 


IIAA'avglh 

) 


\\K cross 1 1 2 




Case 1: f) -- 


= 10" 8 , e = 


10-", M surf = 1 




Set 1 


0.2415 


1.63x10" 


" 4 0.9788 


2.261 


Set 2 


0.2415 


1.61x10" 


" 4 0.9737 


2.254 


Set 3 


0.2415 


1.46x10" 


" 4 0.9753 


2.247 




Case 2 


p = o,e = 


0, M surf = 




Set 1 


0.1932 


0.104 


0.4578 


4.253 


Set 2 


0.3439 


0.433 


0.4551 


6.055 


Set 3 


0.2029 


0.134 


0.4495 


5.224 



Does this mean that mixed modes are therefore not 
a useful diagnostic of core conditions, and that previous 
studies somehow exaggerated their importance? Certainly 
not. Rather, it means that the methods presented here are 
not so good at exploiting the full diagnostic potential of 
such modes. Other methods, such as the one presented in 
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iDeheuvels fc Michell (|2011h . or inversions localised around 
the stellar core, are needed to fully benefit from these 
modes. A possible explanation for this failure can be de- 
scribed as follows. Structural kernels associated with mixed 
modes and g- modes, such as the ones illustrated in the up- 
per panels of Fig. (TSl have large oscillatory behaviour near 
the centre of the star, which does not match the target func- 
tion, T, defined Eq. (flTl) . As a result, their contribution to 
SOLA inversions tends to be reduced, as reflected, to some 
extent, in their inversion coefficients. On the other hand, 
their contribution is not reduced in the KBCD approach 
or the (Ais) scaling law, and this consequently leads to av- 
eraging kernels with highly oscillatory behaviour near the 
centre which can also be a source of error. 

8. Discussion 

In this paper we have investigated different techniques for 
estimating directly the mean density of a star. These in- 
clude the usual scaling law based o n the large frequency 
separation, an approach described in lKieldsen et al.l |2008) 
(the "KBCD approach"), and a SOLA inversion technique 
with a non-linear extension. We have shown how all three 
techniques can be put in a similar form, which allows the 
construction of associated averaging and cross-term kernels, 
useful for characterising what is actually being measured 
and for quantifying the error. 

A comparison of the three techniques, applied to the 
sun and to various test cases, shows that the KBCD ap- 
proach and SOLA inversions yield similar results, whereas 
the scaling law based on the large frequency separation pro- 
duces sub-optimal results. The averaging kernels produced 
by the KBCD approach can come close in quality to those 
produced by a SOLA inversion and are substantially bet- 
ter than the averaging kernels associated with the simple 
scaling law based on (Az/), thereby explaining its relative 
success. Hence, if one only wants to estimate a mean density 
and does not need a full diagnostic on the accuracy of the 
result, the KBCD approach should be sufficient. However, 
for a more detailed analysis, the SOLA method seems ideal 
as it is specifically designed to optimise the averaging ker- 
nel while reducing various sources of error. Furthermore, 
possible improvements to the SOLA method, such as in- 
trod ucing weight functions in the cost function Eq. (fT6|) 
(e.g. iRabello-Soares et al.|[l999h . may help to fine-tune this 
method and should therefore be explored. 

As expected from asymptotic theory, near-surface ef- 
fects are shown to have a stronger impact on the (Az/) 
scaling law than on the two other methods, even when sur- 
face corrections are not included in the SOLA approach. 
However, these two methods require a correct identifica- 
tion of the modes whereas the (Az/) scaling law is much 
more tolerant of mode misidentification. Hence, looking for 
an £ and n assignment which provides the best agreement 
between the simple scaling law and one of the two other 
methods can be a useful way of identifying modes. This 
approach to mode identifi cation, although dif f erent from 
the methods prese nted in iBedding fc Kieldsenl (|2010D and 
IWhite et al.l (|201ll i. is also based on the same principles. 

Quantitatively, a precision around 0.5 % or better was 
reached for the sun, for a Cen B (as based on the spread 
of the results), and for various test cases using SOLA in- 
versions (with P = 10" 8 and 9 = 1CT 2 ) and the KBCD 
approach, whereas the scaling law was less accurate with 



a precision of 1 % or worse. However, when applied to two 
stars observed by CoRoT, the spread in values was larger, 
up to ±2.5 % for the first two methods, and ±3.5 % for the 
(Az/) scaling law. The reason for this remains an open ques- 
tion, but does not seem to stem from the fact that these are 
observed stars, as indicated by further tests involving mod- 
els instead (results not included) and the results obtained 
for a Cen B. Rather, it may be that in the part of the HR 
diagram which contains these stars, there are larger varia- 
tions in the density profile with stellar evoluti on. A ±3.5% 
spread is comparable to the spread obtained in lWhite et al.l 
(|201lD fo r stars with T eS bel ow 6700 K. Above this tem- 
perature, IWhite et al.l (|201lf ) obtained larger differences, 
which in some cases reached 10 %, between the true density 
and that obtained from the (Az/) scaling law. This raises 
the question as to how accurate the KBCD and SOLA ap- 
proaches would be in a similar situation, although one also 
has to bear in mind that this temperature range is less 
relevant for solar-like oscillations. One notable difference 
between the result s described here and those presented in 
Whi te et al.l f|201 lh is that the latter are calibrated to solar 
values rather than nearby reference models, which is likely 
to be an additional source of error, especially for stars with 
a large temperature difference with respect to the sun. 

If combined with a 2 % error bar on measured radii, as 
is expected from GAIA, a 2.5% to 3.5% uncertainty on 
the mean density yields an 8.5 — 9.5 % error bar on the 
mass of the star, where the dominant contribution to the 
error comes from the uncertainty on the radius. Such an 
error bar will then directly translate itself into a compara- 
ble uncertai nty on the inferred m ass of any potential exo- 
planet (e.g. iGuillot k, Have]|[2011l and references therein). 
Accordingly, one may try to obtain more accurate results 
through more complicated approaches such as searching for 
best fitting models in a grid or applying non-linear inversion 
techniques. Such approaches will overcome the non-linear 
effects shown Fig. [5] and which are beyond the control of 
the methods presented here. However, given the difficulties 
in optimising the averaging kernel beyond a certain point 
for a given set of modes, one is still left with the question 
as to how much information is intrinsically available in the 
pulsation frequencies and how much is being assumed in 
these more sophisticated methods. 

Finally, it was shown that the methods presented here 
do not fully exploit the diagnostic potential of mixed modes. 
Indeed, mixed modes have kernels with highly oscillatory 
behaviour near the centre which strongly differs from the 
target function, T, defined in Eq. (|17l) . This reduces the 
contribution of mixed modes to SOLA inversions and leads 
to averaging kernels with an oscillatory behaviour near the 
centre for the (Az/) scaling law and the KBCD approach. 
Also, the linear domain for these modes is smaller and 
can lead to substantial errors when outside this domain, as 
shown in Sect. 16.51 However, one may expect these modes 
to be very useful when applying localised inversion methods 
to models which are close to the true stellar structure. 
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